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O ■ Abstract 

q-i There is proved an existence theorem, in the Newtonian theory, for 

static, self-gravitating, isolated bodies composed of elastic material. The 
bi). theorem covers the case where these bodies are small, but allows them to 
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have arbitrary shape. 
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1 Introduction 

Most solutions of the Einstein field equations — whether known explicitly or 
given by existence theorems — describing static isolated bodies are spherically 
symmetric. The reason for this is the choice of matter model which usually is that 
of a perfect fluid - and such models are necessarily spherically symmetric. This 
latter statement has been proved by Lichtenstein |7| in the Newtonian theory and 
- in the same generality - is still a conjecture in General Relativity. (For the best 
results available, see ||, 0). The only nonspherical solutions known to us are 
the axially symmetric ones constructed by Rein for Vlasov Matter [|ll|]. In the 
present work we pursue another way to describe nonspherical gravitational fields 
by choosing as our matter model elastic bodies, coupled to the static Einstein 
equations. In the spherically symmetric case this has been done by Park [HJ. In 



the nonspherical case nothing is known in the Einstein theory nor — to the best 
of our knowledge — in the Newtonian theory. Thus, as a first step, in the present 
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paper we prove an existence theorem, purely in the Newtonian theory, for static 
self-gravitating bodies composed of elastic material. The theorem allows these 
bodies to have arbitrary shape. 

The main limitation of the present work is that we restrict ourselves to solutions 
close to the natural state of the body, which, in physical terms, means that we 
require the body to be "sufficiently small" . The main technical tool, then, is the 
implicit function theorem near that natural state. For "pure traction problems" 
such as the one studied here, there often occur phenonema of non-uniqueness 
beyond the trivial one stemming from invariance under Euclidean motions. These 
phenomena, which have been thoroughly studied (see Chillingworth et al 0) 
do not happen in the problem at hand. The reason is that traction problems 
require, e.g. in the case of vanishing traction, certain compatibility conditions 
("equilibration conditions") on the load, namely that the total force and total 
torque it exerts on the body be zero. In our case, where the load is given by the 
pull of the body's own gravitational field, these quantities are a priori zero. 
We also point out that we are not able to make statements on the global problem, 
i.e. what happens far away from the natural state. For this one would invoke 
variational techniques, in particular the powerful methods introduced into the 
subject by Ball (see e.g. 0). 

Let B be an open, bounded, connected subset of R 3 with smooth boundary. The 
domain B ("body") is our reference configuration. We also consider 1-1 maps 
: B -> R 3 ("physical space"), x { = 4>i{X A ). Let B^ C R 3 be the image of B 
under (/). Then the basic field equations are as follows: 

- div x T* = p grad x U in B^ (1) 

A X U = AnGp in R 3 (2) 

Here G is the Newton constant, U is the gravitational potential, is the sym- 
metric Cauchy stress tensor, the mass density p satisfies p = np with p a positive 
constant in B* and po = i n R 3 \ ~B ■ Finally, the number density n : B* — > R 3 
is given by n(x) = detV/(x), where / is the inverse map of <fi. We assume 
n > in B*. Let us remark that Equ.([l]) and Equ.@ also describe perfect fluids, 
namely if it is assumed that = pi, where the pressure p is a function just of 
n and Ijj = Sy. 

We first make some observations on the equations ([!]) and @ separately. The 
divergence structure of (|1|) implies that its right-hand side, say b, satisfies a 
compatibility condition, as follows: Let £i(x) be a Killing vector of R 3 , considered 
as flat Euclidean space, i.e. of the form 

£i(x) = Ci + UijXj <^> di€j + dj£i = 0, (3) 

where q and u>ij = u)[ij] are constants. Then, upon scalar multiplication of ([I]) 
with £ and integrating over B*, we easily find that 
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[b>t+ I ^-£ = 0, (4) 

JB* JdB<f> 

where tt = T^uf with the outward unit-normal of dB*. 
On the other hand the "load" b inserted on the r.h.side of ([!]) has the property that 
it gives no contribution to Equ.(|J). This is seen as follows. Define a symmetric 
2-tensor by 

1 1 

@ = -^{gradJJ <g> gradJJ - -I gradJJ ■ gradJJ). (5) 

Then Equ. (0) implies 

div x Q = p gradJJ (6) 
Suppose, in addition, that U satisfies 

U = 0(A) on R 3 \ (7) 

\x\ 

Operating with £ on @ as before on (|l|), but with integration over R 3 , we find 
that the load b used in Equ.(|l]) is "automatically equilibrated" in the above sense. 
Put differently, choosing for £ the three translation Killing vectors, this statement 
amounts to saying that the force exerted on the body by its own gravitational field 
is zero. Similarly, using three rotational Killing vectors, implies the vanishing of 
the gravitational self-torque. 

We want to solve the coupled system (ED) and (0) subject to no-traction boundary 
conditions, namely that tf be zero on dB^, which is a free boundary. To make 
the problem tractable it is thus important to write the above equations as PDE's 
on B, rather than B^, using the Piola transform. With the definition Tia = 
n ~ 1 fA,jTfj, one finds (see e.g. ||) that 

— divxT = pq grad x U. (8) 

If Tfj is solely a function of f^, TiA can be viewed as a function of (V0)^ = F^a- 
This follows from the chain rule for differentiation. The potential U(x), satisfying 
(H) and ([?]), is given in physical space by 

U(x) = -G Po I , H ^ , d 3 x'. (9) 

JB<i> \X — X'\ 

Consequently (§) takes the form 

MX) -MX') j3v , 



- d ^ G *L \; ( x)4(x>w ix ' (10) 



and the compatibility conditions @, using (|3|), result in 
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dB 



t=0 (11) 



tA(/)(x) = 0, (12) 

where U = T^aVa- Our aim is to solve ([LC]) for 0, subject to the boundary 
conditions 



\8B 



(13) 



We assume that 
Al T iA (V(j)) = 0, when V<p = I 

A2 The linearization- at- (0 = id) of the operator divxT is strongly elliptic, in 
other words, a iA jB = d p satisfies a iA jB = a jB iA and aiAjB\F=iViVjVAV B > 


The physical meaning of condition (A2) is as follows: The natural state is usually 
supposed to be such that (iiAjB\F=i = M^ij^AB + Sa^ib) + ^$iA$jB for constants 
/j, and A, the Lame moduli. The ellipticity condition (A2) is then equivalent to 
the inequalities fi > 0, 2// + A > 0. A different interpretation of (A2) would be by 
saying that plane waves propagating according to the linearized-at-F = / time 
dependent equations have real frequency. 

We note that (at least) condition (A2) rules out fluids. And, indeed, the theorem 
of the next section stating the existence of bodies of arbitrary shape, can not 
possibly apply to perfect fluids, as noted in the Introduction. 



2 The Main Theorem 

We now state our precise assumptions. As configuration space C we take maps 
: B -»• R 3 with ^ G W 2 > P (B) 3 , p > 3, and in it C e C C of maps = X { + h { 
with || h \\2,p< e - For e sufficiently small, <fi is C 1 - map close to the identity with 
C 1 - inverse (see Appendix). For the stress tensor Tj^(V0) we assume that it is 
in C 2 (R 9 ,R 9 ) and that it satifies conditions (Al,2) of section [3], wherefrom it 
follows [13] that the operator i— > T(V0) is a C l - mapping from W 2,P (B) 3 to 
W /1,P (B) 9 . Our main result is 

Theorem: For sufficiently small G there is a solution <p G C e of ( |I0D subject to 
(O). This solution is unique provided 



hi{0) = 0, ^(0) = 0. (14) 

(We assume that G B.) 
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We remark that the smallness-condition for G can of course, by scaling, be 
rephrased by Gp <C where L is a typical length scale of B and |T| an 
upper bound for the stress tensor. Our method of proof follows the geometrical 



treatment of the Stoppelli theorem |12[ due to LeDret ||. 

Proof: Consider the map E : C h+ y = {(&,*) G W P (B) 3 x W x - x/pjp {dBf,p > 
3}, the operator of nonlinear elasticity given by 

^ {-div x T{V(j)),T{V(j})v) (15) 

The operator E is well-defined and C 1 (see ||). Now recall from the discussion 
of Section [I] that elements (b, t) G y lying in E(<p) satisfy the compatibility 
("equilibration") conditions, namely 

b+ [ t = (16) 

B JdB 



bA(j)(X) + / tA6(X) = (17) 

JdB 

The set 3^> of pairs (b, t) G y satisfying ([TBI) and (|l7f ), for given 0, is a vector 
subspace of y of codimension 6, when G C e . (More precisely, it follows from 
the results of LeDret, and is easy to check, that the only elements G C, at 
which y^ fails to have codimension 6, are those for which the image 0(B) is 
parallel to a fixed direction v G R 3 - which is impossible if <p G C e .) Let us 
choose some complement S of L e C y e , where L e = y id and define a projection 
P '■ y^> l— ► L e . The linear maps P : y^ i— > L e are isomorphisms and C 1 (see 
P, proof of Proposition 1.4). Next consider the ("live") load afforded by the 
gravitational force, i.e. b = GU((f>) with GUi(4>) given by the right-hand side of 
(|10D. By explicit calculation, or from the discussion of section]!], it follows that 
U ((f)) C 3V Note that this requires B to be connected. If B had several connected 
components, U(4>) would be automatically equilibrated only with respect to the 
whole of B, whereas flT5| ) and ( |T7D would for the operator E be required to hold 
separately for each connected component of B. It is thus important that we have 
only one body. 

We want to solve the equation 

E((j)) = GU(4>) on B, U = on dB (18) 
for small G. We know from Al that = id is a solution for G — 0. We write 

F(G,4>)=E(cf>)-GU( ( t>), (19) 

with F viewed as a function R x C ( h J. In the Appendix we show that U, 
whence F, maps C e into 3^ in a C 1 - fashion. If we now compute the linearization 
of F at = id for G = 0, we find that this is a map from C to y which is not 
surjective, due to the presence of the equilibration conditions. To get round this 
difficulty, consider the modified operator 
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F'(G,<t>) = P(E( ( f>))-GP(U( ( f>), (20) 

with F' viewed as a map R x C e i— ► L e . Clearly every solution of F' — is also a 
solution of F — 0. If, in addition, we eliminate the translational and rotational 
freedom by replacing C by C sym , consisting of all elements 4>i — Xi + hi in C for 
which Ui(0) = and duUjAO) = 0, it follows from standard linear theory (see 
0, Lemma 3.17 of Chap. 7), that the linearization- at- (0 = id) of F' at G = is 
an isomorphism C S3/m i— > L e . Hence our claim follows from the implicit function 
theorem. 



3 Appendix 

The Poisson integral used in the body of the paper is given by 

Here B C R 3 is bounded, not empty, open and connected with <9B smooth and 
4>i G W 2 ' P (B) 3 , p > 3. Furthermore 4>i(X) = Xi + hi(X), and we assume that 
II ^ 1 1 2,jp< e> e small. It follows from Sobolev embedding that djhi is small , in 
particular bounded in B. By the mean-value theorem we infer that 

\h{X) - h{X')\ < C'\X - X'\ (22) 

This is immediate for B convex, otherwise see H, p. 224. Making e smaller, if 
necessary, we have that 

\dh{X)\< 1 - (23) 
Consequently, there exist positive constants E, E' such that 

E\X - X'\ < \4>{X) - (j)(X')\ < E'\X - X'\ (24) 

It immediately follows that Ui[<f>] G C°(B) 3 . Thus Ui is a bounded map of 
U e (id B ) C W 2 ' P (B) 3 -> C°(B) 3 , whence to to W°' P (B) 3 . Here U e (id B ) denotes 
the set of 0's in the e-ball centered at the identity map. We want to show Ui is 
actually C l . We first compute the Gateaux-derivative (directional derivative) of 
Ui, namely 



*W1 



DUi[cjP\-v (25) 



t=o 



where 4>\{X) = Xi + G W 2 ' P (B) 3 . First observe that the expression 
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1 f MX) - g(X) - 



(26) 



for 1^1', converges pointwise for t — > to 



- «,(X0 3(#(JT) - ^(X')M(X) - 0O(X'))(^(X) - Vj (X')) 



<P(X) - <p°(X>)\ 3 



p(X) -0°(X')| 5 



(27) 



Next note the following chain of elementary inequalities: &i,&o vectors G R 3 (or 



1 1 K I61-60I 



1 1 , / 1 1 \ „ , . 



N 2 N 



N 2 N \b \ 2 \h\ 



' N 3 N 3 
ai, 6i, oo, ai vectors G R n 



1 1 , / 1 1 1 \ „ 

< I + — ^rr, ~ + t^tt^ I \h - 6o| 



N 3 \h\ 



< 



N 3 N \bi\ 2 \b \ 2 \h\\b \ 



1 1 

+ 



|6x| 3 |6 | N 2 M 2 NN 3 , 
\ai - a \ 



+ 



M 3 



setting 



a l = b l = <p t (X)-<f) t (X > ) 
a = b = <f>°(X) - 4>°{X') 



we find 



(28) 
(29) 
(30) 



|ai||6i-6d| (31) 
(32) 



(33) 
(34) 



tip) - #m Aw - Am 



|^(X)-0<(X')| 3 |0°(X)-0O(X')| 3 
( 2 

+ 



< 



+ 



b°(x) -0°(x')| 3 \<p(x) - <p(X')\ 2 \<iP(X') -(jP{x)\ 
i 



\Mx>) - <pt{x)\\nx>) - nx)\\ 



t \v(X')-v(X)\ 



(35) 



It follows, using (p4|), that the sequence in (p6| ) is bounded by a positive, £- 
independent function, whose integral over X' G B is a bounded function of X 
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in B. So, by dominated convergence, the previous limit, whence the Gateaux 
derivative actually exists. But the linear operator defined by the directional 
derivative v € W 2,P (B) 3 — > W 0,P (B) 3 , is clearly bounded. So, by a standard 
theorem (see e.g. Corollary 2.4.10), U(<f>) is a ^-functional. 

Acknowledgments: The authors thank Horst Beyer and Alan Rendall for useful 
discussions. 
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